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Abstract. The motion of circular WP for one electron in central Coulomb field 
' with high Z is calculated. The WP is defined in terms of solutions of the Dirac 

equation in order to take into account all possible relevant effects in particular the 
^^^^ ^ spin-orbit potential. A time scale is defined within which spin dynamics must be taken 

^ ' into account mainly in the atoms with high Z. Within this time scale there exists 

a mechanism of collapses and revivals of the spin already shown by the authors for 
' harmonic oscillator potential and called the spin orbit pendulum. However this effect 

. has not the exact periodicity of the simpler model, but the WP's spatial motion is 

' nevertheless quite similar. 

o 
o 

:^ 
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qh| 1. Introduction 



In the last decade large efforts |jT]-0 have been devoted to the detailed understanding 
of the quantum dynamics of WP in simple systems like H atoms or hydrogenoid 
atoms. Theoretical investigations resulted in very good understanding of such subtle 
interference effects like collapse, revivals, fractional revivals and super-revivals of 
wavepackets created in atomic and molecular systems. Many signatures of these 
phenomena have been observed in numerous sophisticated experimental studies. Indeed 
with contemporary lasers and tunable electromagnetic fields one can engineer many 
different desired initial states and analyze their time evolution ||^. 

Our aim in this paper is to perform a new step in the analysis of the WP dynamics 
for hydrogen atom and for ions with one electron by considering spin motion. Until now 
this effect was not considered for isolated atoms since it requires very long time scales 
to manifest itself due to the weakness of the spin-orbit coupling. In order to decrease as 
much as possible this time scale we will consider WP within ions with very large Z. We 
intended to make a calculation which include all possible effects for these elements and 
therefore we have considered WP which evolve under the relativistic Dirac Hamiltonian 
i.e. which are four component spinors. 
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Spin effects and relativistic effects have been studied many times but mostly for 
atoms submitted to intense laser fields. A recent review has been given in [|1^]. Some 
of the most relevant research works are found in |13-16|. It is generally agreed in them 
that spin and relativistic effects are small but that their magnitude increase with the 
intensity of the field. For fields of the order of 10^^ to 10^^ W/cm^ they should be taken 
into account. In a most recent paper [W\ the interest has been focussed on the spin- 
orbit coupling, essentially for Al+^^ and Ga+^^ ions. The dynamics of the spin degree of 
freedom was analyzed there during the time of interaction of the laser with the ion. The 
authors did not use the Dirac equation but the Schrodinger equation corrected up to 
order and therefore the spin-orbit potential could be clearly identified by its mere 
suppression. The most significant signature occur in the radiation spectrum. However 
little polarization was produced for the ions which were considered. 

Spin dynamics was previously studied by us in harmonic oscillator potentials in 
non-relativistic dynamics. In ref. [p|-[TT[l we have found a new effect that we have 



called the spin-orbit pendulum for a spin-orbit potential with a constant radial part. 
This phenomenon is simply the periodic collapse and revival of the average spin and 
an exchange of angular momentum between spin and the orbital angular momentum. 
The analysis is quite parallel to that of the Jaynes-Cummings model In order to 
extend this situation to a relativistic dynamics it was quite natural to study the Dirac 
oscillator (DO). There, the strong entanglement that gives rise to spin oscillation is still 
present [0. Because of nonlinearities of the relativistic energies the dynamics is not 
periodic. We have also studied relativistic effects like zitterbewegung. Our studies have 
been also extended to WP motion corresponding to linear trajectories in non-relativistic 
approach (HO+LS Hamiltonian) [|1^] as well as in the relativistic Dirac oscillator [pO|] . 

Our previous studies stimulated the present work for hydrogen or for heavy ions. It 
is expected that many features of the DO should take place in this much more realistic 
situation. 

Our paper is organized as follows. In order to make it self-contained we in Sec. ^ 
present basic formula for solutions of Dirac equation with Coulomb potential. In 
Sec. D we show the construction of the circular WP and in Sec. | we discuss its time 
evolution, focusing on the autocorrelation function, spin expectation values and spatial 
motion. Particular attention is paid to relativistic effects and different time scales. The 
concluding remarks are contained in Sec. |^. 



2. Dirac equation with central potential 

Solution of Dirac equation for a central field, in particular for Coulomb potential is one 
of standard, textbook problems in quantum mechanics [pT| , |22[] . Therefore we will recall 
here only few formula that will be relevant for our further considerations. 
Dirac equation for a spherical potential V{r) is 

HD^{r) = (ca ■ p + (3mc^ + 1/(r))^(r) = E^(r) . (1) 
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Its eigenfunctions can be found in a representation that diagonalizes simultaneously 
with Hf) three other operators, j'^,jz and K = h + L ■ cr. Then we may write for \&(r) 
(using phase convention as in 



with 



I' = 2j - / 



g{r)^ijm 



2(/ + !)-/ = / + ! for J = l + \ 



2(1-1) -1 = 1 



for j = I 



3+ 



The explicit forms of the spherical spinors ^ijm are: 

I \l~ 



l,j+,m 



j-tn Y 
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and 
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Radial wave functions are given in terms of confluent hypergeometric functions F{a, c, x 



, a a(a + 1) x 



Explicit formula are [^]: 

JXr) = 
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X [n'F{-n' + 1, 27 + 1, 2Ar) - {k - ^/X)F{-n, 27 + 1, 2Ar)] 



and 
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r(27 + 1) 



r(27 + n' + l)(l + E) 
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X [-n'F{-n' + 1, 27 + 1, 2Ar) - (/« - ^/X)F{-n', 27 + 1, 2Ar)] , 

where ^ = 7 = V^^^, = (j + l/2)^ = [1 + and A = v/r^:^^, 

a = /he = 1/137.036 being the fine structure constant. The nonnegative integer 
number n' is related to n and |k| by n = n' + In the following calculations we use 



atomic units 



h = c 



3. Construction of wavepacket 

In our first paper on relativistic WP in Dirac oscillator (DO) we were considering 
a circular spinor without small components written as 

/ a\m) \ f a\ 
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(7) 
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Here a and b are coefficients which determine the initial direction of the spin. There is 
no loss of generality to choose them real. i?„=;(r) denotes radial wave function of the 
harmonic oscillator. The large component for the partial wave / is written as 




3+,3+ 



+ b 



(8) 



Weights Wi of the superposition (|^) were chosen in a way providing \& to be the HO 
coherent state. For DO the bispinors without small components were constructed as a 
superposition of spinors with positive and negative energies that provide a cancellation 
of the small components. 

For the Dirac equation with Coulomb potential we are not able to get rid of small 
components completely, as we want to build the WP only from bound states. We can, 
however, construct small components in a way similar to the big ones providing their 
circular motion. For this purpose we write the angular parts of the large components 
as in eq. (Q), then we complete each of the spinors and by their 

corresponding small components associating each with its radial part different for the 
states with and j_ as in eq. (^. Thus there are four radial parts respectively called 
gj+, fj+,gj- and /j_ for a given I. In order to produce circular waves for hydrogenoid 
atoms the connection between n and / is / = n — 1. 

A partial wave at time t can be written with these notations using Ej_^_ and Ej_ the 
eigenvalues of the Dirac equation: 



+ 



3+ "'+i,i+J- 



-lE+t 



+ b\ 
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ifl'i„ ^i,j-,j- 



-lE-t 



(9) 



After expansion of the Vl as in 
written as 



the circular wavepacket at time t can be 



/ 



(10) 



\ciit)) \ 

\c2it)) 

\c-m 
V ic4(t)) ; 

For Coulomb potential weights w„ (that can be made real) are chosen in a standard 
way 0,^|2^, i.e. as a Gaussian distribution with respect to the mean value A^: 
|^i;„|2 = (27r4)-5e-("-^)'/24_ 

Explicit form of the wavepacket (|1^) is obtained in the following form: 



(11) 
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\c2{t)) = wib\l,l) ^9j+^f^ 
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For t = spin and orbital angular momentum are only approximately decoupled, 
they were exactly decoupled in (|^). 

In numerical calculations we use the full bispinor (12q- 124 ). However, particularly 
for > 10, the contributions arising from the small components {{c^) and {c^}) are 
almost negligible, because the radial /(r) functions are always few orders of magnitude 
smaller than the corresponding g{r) functions. This is due to the fact that their 
prefactors differ mainly by the factor Vl — E (few orders of magnitude smaller than 
1) for /(r) and \/l + E (close to ^/2) for g{r). 

Using eqs. ( |12cj - [L2dD one can calculate the contributions from small components of 
our circular WP exactly. One obtains 

i^^^)smaU = (C3(t)|c3(t)) + {c^{t)\c^{t)) (13) 

where 



(c3(t)|c3(t))=E 
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+ 6^ 



(c4(t)|c4(t))=E^? 
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2/ + 3 



+ b' 
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In fact, {'9^'^)smaii (C3IC3), and (C4IC4) are constant contributions for our WP. Terms Fj_ 
and Fj_^ are radial integrals defined in next section. Fig. |l| displays the magnitude of 
{'^^'^)smau as function of Z and for Z G [1, 92] and G [2, 60]. Only for very large 
Z and very small this contribution exceeds 1%. 



4. Wavepacket evolution 

4.1. Time scales 

We can define several characteristic times on the basis of the previous works about 
recurrent wave packets ■ Writing the energy in function of a single quantum number 
n we define for n = N a. hierarchy of times 
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Figure 1. The magnitude of the smaU components contribution as function of Z 
and N. 



For k = 1 we obtain the classical Kepler time, for k = 2 the revival time, for k = 3 the 
superrevival time etc. .. For hydrogenoid atoms those terms increase as N^^"^. These 
times can be compared to a characteristic spin orbit time Tis = 271 /uis with uis defined 
by 



hujj 



E. 



■j-=N -3/2 



(16) 



:i7) 



To first order we obtain 

^ 1^ ■ 

Therefore the spin-orbit frequency follows a variation parallel to the superrevival time 
with k = 3. Those times are plotted on Fig. H for N > 2 and for a large range of 
atoms (ions) including uranium. For hydrogen Tig is orders of magnitude larger than all 
the other characteristic times T(l) — T(5). If we want to observe spin-orbit effects in 
the picosecond regime we can concentrate on = 20 for Z = 92. With such a choice 
Tu ^ 1685 T(l). The other times are T2 = 13.33 T(l), T3 = 200 T(l), T4 = 3200 T(l). 
The choice of Tig as a convenient time is justified by our plot of the autocorrelation 
function of Fig. ^ which shows an interesting oscillations up to 12 periods Tis. 



4.2. Autocorrelation function 

The autocorrelation function (often referred as the recurrence probability) contain a lot 
of information on wavepacket evolution. It is a convenient measure of the degree of 
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recurrences and can be defined as: 



Ait) = i^m^it)) = J:{c^mcm . 

i=l 

Using equations ( |12a| - P-2c|) one obtains 
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+ 6^ 



2/ + 3 



F,-+e-'^^+* 
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are radial integrals (for instance Gj^ = r"^ g^^dr 



and Gj^ = /q°° r gj^gj_dr and so on) tliat can be easily expressed in terms of F functions. 

In the presentation of our results below we focus (as in our previous papers) on 
the most interesting case, i.e. a = b = l/y/2. This means that initial direction of the 
spin is along the Ox axis, in the orbit's plane. In Fig. ^ we present the square of the 
autocorrelation function |A(t)P as function of time for long term evolution (T;^ units) 
of the WP with = 20. Heavy ion system with Z = 92 is considered with the full 
relativistic approach. Very regular revival structure connected with integer values of 
time in T;^ units is immediately recognized. This regularity is somewhat distorted by 
a strong influence of T4 periodicity (see right bottom picture of Fig. ||) but becomes 
even more evident for larger N. Both the main part and the insert of Fig. ^ show very 
complicated shorter scale behaviour of the autocorrelation function related to the three 
shorter time scales: Tci,Trev,Tsr, known already from the non-relativistic studies 0. 
The case with = 40 shown in Fig. ^ suggests that for times t G [0, 1.5]T;s and 
t G [6, 10] Tis the T;^ time scale is essential for the evolution of the WP. It shows also 
that there exist a deep connection between the time behaviour of the autocorrelation 
function and the average values of the spin operators as discussed in the next subsection. 
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4-3. Expectation values of spin operators 

Expectation values of spin operators are derived directly from ( |12a| - [12Q| ). One obtains: 



'y 5 



(205) 



{2/ 1 QJ 
^'G,, + b'j^^.iG,, - 21G,_) - h'j^^^.G,, cosu.t 

+b -. F. -\a^—^ + b^— — F,- ), 20c 

21 + 1 ^' \ 21 + ?, (2/ + l)(2/ + 3); ' ^ ' 

where uji = {Ej_^ — Ej_)/h. Terms denoted by Sa^^ and 6ay, which are only small 
contributions, stem from coupling between states with / or j different by 2 units. Their 
explicit expressions are: 

da^ = '^Ki cos LJit , (21a) 



where 



6ay = ^^Ki sin uit , (215) 
I 

huJi = Ej^ - F(j+2)_ , (22 a) 

F'= / f)+fa+2).dr , (225) 



and 



Ki = 2abwiWi+2F\ 



(2/ + 2)(2/ + 4) 



\\ (2/ + 3)(2/ + 5) ■ 

In a non-relativistic limit the terms containing the radial parts fj^ and fj_ cancel 
and the correction terms Sa^ and 6(jy tend to zero, the same for Fj_^ and Fj_. For the 
harmonic oscillator with a constant spin-orbit potential the radial integrals verify 

= G,_ = 1 . (23) 

The dynamics of the spin motion in such a case has been discussed by us in and |[T0(1 . 
Using the dispersion law appropriate to a constant spin-orbit potential 

uji = {21 + l)ui, (24) 

we have obtained there the periodic oscillations between spin and orbital angular 
momentum that we have called spin-orbit pendulum. 

In the present, more general case we see that both (a^) and (a^) contain an 
additional constant shift with respect to the non-relativistic expressions, coming clearly 
from the small components. In addition there are time dependent shifts due to the 5ax 
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and 6ay terms. Numerically these shifts may be safely neglected. Therefore the most 
interesting dependence comes from the use of the eigenvalues of the Dirac hamiltonian 
in the time behaviour that we will now discuss in more details. 

In Figs. H and |^ we present the average values of the spin operators and the length 
of the spin vector as functions of time for cases N = 20 and N = 40 [Z = 92). The 
latter case the most distinctly reminds both non-relativistic spin-orbit pendulum for the 



HO+LS Hamiltonian ||y,|IO[ and the relativistic one for the DO |18l]. We clearly recognize 
time ranges for the spin collapse (0 < t < 2.5), strong spin entanglement (2.5 < t < 4.5) 
and spin revivals (6 < t < 10, all times in Tig units). For the case N = 20, shown in 
Fig. we see a similar behaviour of the spin averages, distorted, however, much more 
due to stronger interferences with T4 characteristic time. In cases discussed above the 
ratio T4/T/S is approximately equal to 1.9 for = 20 and 3.7 for = 40. 

4-4- Spatial motion of circular WP 

Let us now discuss the motion of WP in space. A complex structure of the 
autocorrelation function and the existence of several time scales suggest that 
qualitatively similar spatial evolution can occur in different time scales. Below we focus 
our presentation of the probability density motion in two particular time scales: the 
shortest one, corresponding to the Kepler period of classical electron motion on the 
circular orbit with n = N, and the spin-orbit motion time scale. The admixture of the 
components |/, /—I) in |ci(t)) (|12q|) [as shown earlier the small components give negligible 



contributions to causes deviations of the WP motion from the equatorial plane. 



known already from our previous non-relativistic studies |TT| , p^ . The magnitude of these 
deviations decreases, however, when increases and for > 10 — 12 the maximum of 
the probability density moves very close to the equatorial plane, i.e. the plane of the 
classical orbit. Therefore we limit our presentation of {t) to this plane only. For the 
short time scale, t E [0, 9/4] T^, it is shown in Fig. 0. The sequence of pictures exhibits 



the behaviour known from non-relativistic studies on WP in hydrogen atom |]l]-^|^,^ 



namely circular motion with spreading along the orbit and forming interference patterns 
that lead to fractional revivals, n Fig. ^we display the full probability density \I'^\I^, not 
the spin up and spin down components separately, because for such short time scale no 
difference in their motion is visible. 

In order to see this difference one needs to observe WP in time scale comparable to 
Tig. In Fig. I we present the sequence of shapes of the same WP (A^ = 20, Z = 92) for 
times t = 0, 1/2, 1, 3/2 and 2 (in Tis units). In these cases the presentation is divided into 
the two parts: in the left column spin up components of are displayed, whereas in 
the right column the spin down components are shown. At the presented time instants, 
that are multiples of T/s/2 both sub-packets are quite well revived (localized) and move 
on a circle with different velocities. 

Fig. P shows (in the same convention) the shapes of WP's components with spin 
up (left) and spin down (right) for time t = 10.063545 T/^ (A^ = 20, Z = 92). As seen 
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from Fig. ^ this time corresponds to the highest value of the autocorrelation function 
(|y4(t)p ^ 0.8) in the presented range of times. One sees from this figure that the degree 
of the recurrences (of both sub-packets separately) is indeed very high in this case. 



5. Conclusions 



We have discussed the full 3+1-dimensional evolution of the relativistic wave packets in 
Coulomb potential. The main relativistic effect is the appearance of the new time scale 
due to the spin-orbit coupling. In general this effect has the same origin as spin-orbit 
pendulum, introduced by us both in non-relativistic |]9HlI| and relativistic |jl8|] studies, 
where harmonic oscillator potential with a strong spin-orbit coupling was considered. 
In the Coulomb potential (relativistic) case presented here the time scale of the spin- 
orbit motion is large, therefore there is practically no chance for observing it in such 
systems like hydrogen or hydrogenoid atoms. However, this time scale decreases rapidly 
as function l/Z"^, when heavy ions with the single electron are considered as systems 
for experimental studies. At Z = 82, and even better Z = 92 the Tis time becomes 
reasonably short and in principle experimental observations of the spin-orbit effects in 
electron WP's motion may become possible with existing techniques. This conclusion is 
based on the fact that the single electron heavy ions with Z = 82, 92 have been already 



created and some of their properties have been measured [|25|. It is also possible to 



suitably adjust the two parameters N, Z in order to fit the energies of electron states 
and time scales of the motion to ranges of action of available lasers that have to be used 
for creation of wave packets and analysis of their motion. 

Comparing the phenomenon of the spin-orbit pendulum discussed in our previous 
studies jgHTTi |l8|] with the present Coulomb potential case one sees the following 
differences. During their motion Coulomb wave packets experience several important 
time scales (only one is present in the HO case) that make the motion much more 
complex and obscure the spin-orbit motion. Relativistic (nonlinear) dependence of the 
eigenenergies on quantum numbers destroy exact periodicity of the motion. 
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Figure 2. Characteristic times of the WP evolution as fmictions (in log log plot) 
of the mean value of principal quantum number for different values of nuclear charge 
(from Z = 1 to 92). The period of spin-orbit motion is given by the thick solid line. 



Spin dynamics of WP in strong Coulomb field 13 

1 1 1 1 1 1 1 1 1 1 1 1 r 



0.8 




Figure 3. Square of the autocorrelation function for system with N = 20, Z = 92 for 
long term evolution. In this case Tig « 1685 T^/, where = 2tt /{ZaY is a Kepler 
time for non-relativistic WP with the same N . Short term behaviour around the best 
revival is presented in the insert. 
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Figure 4. Square of the autocorrelation function for system with N = AQ, Z = 92 for 
long term evolution. In this case Tis ~ 6921 Td- 





Figure 6. Time evolution of the average values of spin components (cTj), i = x,y,z 
for system with N = 40, Z = 92. 
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Figure 7. Spatial motion of the WP with N = 20, Z = 92 in the initial stage of 
time evolution t G [0, 9/4]Tc;. The probability density is plotted as function of 
the position on the plane of the electron's classical orbit. The radius of the thick circle 
corresponds to the initial distance of the WP's peak from the center of the potential. 
Note that in this figure and all next plots the contour lines correspond to the same 
values of the probability density. 
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Figure 8. Shapes of the WP {N — 20, Z — 92) components with spin up (left 
column) and spin down (right column) during evolution in time scale measured in Tis 
units {t e [0, 2]Tis). The different angular velocities of the motion of these sub-packets 
are clearly visible. 



Spin dynamics of WP in strong Coulomb field 



17 




Figure 9. The best revival of the WP for the system with = 20, Z — 92 presented 
on the plane of the classical orbit. The time instant t = 10.063545 T;^ corresponds to 
the highest value of the autocorrelation function displayed in Fig. ^. For comparison 
with the spin components at t = see the upper row of the Fig. ^ 



